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Abstract 

We prove the existence of Kolmogorov-Petrovsky-Piskunov (KPP) 
type traveling fronts in space-time periodic and mean zero incompress- 
ible advection, and establish a variational (minimization) formula for 
the minimal speeds. We approach the existence by considering limit 
of a sequence of front solutions to a regularized traveling front equa- 
tion where the nonlinearity is combustion type with ignition cut-off. 
The limiting front equation is degenerate parabolic and does not per- 
mit strong solutions, however, the necessary compactness follows from 
monotonicity of fronts and degenerate regularity. We apply a dynamic 
argument to justify that the constructed KPP traveling fronts prop- 
agate at minimal speeds, and derive the speed variational formula. 
The dynamic method avoids the degeneracy in traveling front equa- 
tions, and utilizes the parabolic maximum principle of the governing 
reaction-diffusion-advection equation. The dynamic method does not 
rely on existence of traveling fronts. 



*Department of Mathematics, University of Texas at Austin, Austin, TX 78712 
(jnolen@math.utexas.edu). 

^Department of Mathematics, University of Texas at Austin, Austin, TX 78712 
(rudd@math.utexas.edu). 

* Department of Mathematics and ICES (Institute of Computational Engineering and 
Sciences), University of Texas at Austin, Austin, TX 78712 (jxin@math.utexas.edu). 



1 



1 Introduction 



Reaction-diffusion front propagation in heterogeneous fluid flows has been 
an active research area for decades ([I], [T3], [T3], |22|, [23- [H2] and 
references therein). A fascinating phenomenon is that the large time (large 
scale) front speed is a highly nontrivial quantity that depends on the multiple 
scale structures of the fluid flows. Speed characterizations have been studied 
mathematically for various flow patterns by analysis of the reaction-diffusion- 
advection (RDA) equations (see [1131311111111113113111112011231231^1211 
and references therein). In particular, variational characterizations have led 
to speed asymptotics in both deterministic and random flows (23 H3 E3 123 

I22H23- 

In this paper, we consider front solutions to a RDA equation of the form: 

u t = Am + b(x,t) ■ Vw + f(u) , (1.1) 

where u = u(x, t), x G ~K N , and t G R. The N components of the vector field 
b(x, t) := (ft 1 (a;, t), b 2 (x, t), . . . , b N (x, £)) are smooth and spatially divergence- 
free, are periodic of period 1 in both x and t, and have mean zero over the 
period cell Q x (0, 1), where Q is the unit cube in R^: 

V x -b = and / W(x, t) dx dt = 0, for i = l,...,JV. (1.2) 

JQx(0,1) 

The function / : [0, 1] — > R is KPP, namely a C 2 function of u such that: 

/(0) = f(l) = 0, /'(0) > 0, /'(l) < 0; /(«) > , Wu G (0, 1); /(«) < w/'(0). 

(1.3) 

An example is f(u) = u(l — u). 

A more general class of nonnegative reaction nonlinearities arising in ap- 
plications may allow degeneracy at zero, e.g. f(u) = u m (l — u), for m > 2, 
and f(u) = e~ E l u {\ — u),E > 0. Let us define the so called positive non- 
linearity to include KPP and generalizations. A positive nonlinearity f(u) 
satisfies: /(0) = /(l) = 0; f(u) > 0, u G (0, 1). In the course of this study 
we will also consider the ignition type combustion nonlinearity, for which 
f(l) = 0, /'(I) < 0' an d there exists a constant 6* G (0, 1) such that 

/(u)=0 for ue[Q,0] and f(u) > for ue(9,l). (1.4) 

The constant # is known as the ignition temperature. 
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Given a unit vector k = (k 1 , . . . , k ) G K , we look for a planar traveling 
front which propagates in the direction k with constant speed and solves 
(jl.ip . Due to the form of the advection field, we seek a solution of the form 

u(x, t) — U(k ■ x — ct,x,t) — U (s, y, r) (1.5) 

for s := k ■ x — ct, y := x, and r := t. The constant c is the speed of the 
front. We require U(s,y,r) to be 1-periodic in both y and r and to satisfy 
the following conditions at s = ±00, uniformly in y and r: 

lim U(s, y, r) — 1 and lim U(s, y, r) = . (1.6) 

s— »— 00 s— »oo 

Substituting [/ into equation (jl.lj) and imposing these auxiliary conditions 
yields the problem 

' u T = {kd s + v^) 2 u + b ■ (kd s + v„) c/ + cu a + f(u) , 

< U (s, r) 1-periodic in y and r, (1.7) 

lim U(s, y, r) — and lim {7(s, j/, r) = 1 . 

Equation ()1.7|) contains only directional derivatives that do not span V Sj2/)T Z7, 
hence there is a lack of compactness. One can view (|1.7J) as a degener- 
ate parabolic equation, parabolic along (s,y,r) = (— c, 0,1), and degener- 
ate elliptic in (s,y,r) = (k t ,e t ,0), e l the standard i-th unit vector in R N , 
i = 1, 2, ■ ■ • , N. To proceed with existence, we shall regularize (J1.7)) then 
pass to the limit to prove a certain weak solution of ()1-7|) defined below. 

Definition 1.1. A front traveling at speed c / with reaction f is a lo- 
cally integrable function U(s,y,r) G L] 0C (1H x Q x (0,1)) whose directional 
derivatives 

U T -cU s , tfUs + Uyi, i = l,...,N, and (kd s + W y ) 2 U 
are continuous and satisfy the traveling front equation Ji.7| ) in the form: 

U T - cU s = (kd s + V,) 2 U + b ■ (kd s + V„) U + f{U) , (1.8) 
and the associated boundary conditions. 
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Our first main result is that such a front solution exists for positive non- 
linearity. The proof is based on elliptically regularizing (|1.8|) and approxi- 
mating positive nonlinearity by a sequence of combustion nonlinearities with 
ignition cut-off. The regularized problem can be solved by existing methods 
in the literature, and one key property is that U is monotone in s. The 
monotonicity, the analysis of asymptotics near s infinities, and the available 
derivative estimates of (jl.8j) imply enough compactness for passage of limit 
and constructing a desired solution. 

A few remarks are in order on the traveling fronts. In the original co- 
ordinates (x, t) G M. N x M, we see that a traveling front U corresponds to 
the function u(x,t) = U(k ■ x — ct,x,t). This is a classical solution of ([1.1)1 . 
because derivatives in the original (x, t) coordinates translate to directional 
derivatives in traveling front variables (s, y, r). Note also that this definition 
of the traveling front is different from the definition of pulsed traveling fronts 
in time-independent advection fields [2]. In that context, a pulsed traveling 
front in direction k G S N ~ X is defined by the relation 

u(x + e,t) =u(x,t-^—^), V(x, t) G R N+1 , 

c 

for any vector e G R such that b(x + me) = b{x) for some integer m. If 
this relation were imposed on the solutions considered here, we would deduce 
that 

u(x + e,t) = U{k ■ x + k ■ e — ct,x + e,t) 

= U(k ■ x — c(t -),x,t) 

c 

would be equal to 

U{X, t ) = U (k ■ X — C{t ), X, t ). 

c c c 

Hence, if b is T-periodic in time variable t, we would conclude that 

k ■ e 

c = — — , tor some integer mf 
mT 

This implies that if b has a nontrivial time- dependence, the speed c would 
be determined only by the period of b in space and time, regardless of other 
information on b, which is nearly impossible. The price of seeking a more 
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general form of traveling fronts (jl.5|) is the difficulty caused by degeneracy 
in front equation ([1.8)1 . 

Degeneracy makes it hard to prove the minimality of front speed, if one 
relies entirely on the front equation (jl.8|) . We adopt the alternative approach 
of studying the time dependent front-like solutions of (jl.lj) where there is 
dynamics but no degeneracy. The KPP nonlinearity is critically used in 
constructing comparison functions and passing to the large time limits of 
solutions. 

Our second main result is that the KPP traveling front speed is minimal 
(denoted by c*) and is given by the variational principle: 

c* = inf /i(A)/A, 
A>o 

where /z(A) is the principal eigenvalue of the periodic-parabolic operator [T2] 

L A $ := A a .$ + (b - 2Xk) ■ V x $ + (A 2 - Xb ■ k + f (0)) $ - $ f , 

defined on spatially-temporally periodic functions The variational 

principle reveals that the flow structures in b are upscaled into the front speed 
c* via the principal eigenvalue /i, see jH] for a similar formula. 

The rest of the paper is organized as follows. In section 2, we study 
regularized traveling fronts under combustion nonlinearity with cut-off, and 
obtain necessary estimates to remove the regularization, and justify boundary 
conditions at s infinities. Front existence for positive nonlinearities follows. 
In section 3, we study propagation of KPP fronts, minimality of front speeds 
and their variational characterization. In section 4, we present properties of 
the associated eigenvalue problem. Concluding remarks are in section 5. 

2 Existence of Fronts: Positive Nonlinearity 

In this section, we study the problem of existence of planar traveling front 
solutions of ([Lip for general positive nonlinearities. We prove the following: 

Theorem 2.1. Let f be a positive nonlinearity: f(u) > 0, u G (0, 1); f(0) = 
f(l) = 0. Then there exists a traveling front U* with speed c* > that solves 
il.fy) with the boundary conditions in \1. 7| ). 

Our proof of Theorem 12.11 proceeds by first solving a regularized version of 
(|1.7|) where we approximate / by a combustion-type nonlinearity f e such 
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that f — > f as 9 — > 0. In Section l2~Tl we prove existence of solutions to the 
regularized problem for the combustion nonlinearity f 9 . Then, in Section 12 .21 
we show that as 9 — > 0, f e — > / and the solutions converge to solutions of the 
regularized problem for the general positive nonlinearity. Finally, in Section 
12.31 we complete the proof of Theorem 12.11 by removing the regularization 
parameter to recover a desired solution. 

2.1 Regularized Equation: Combustion Nonlinearity 

In this section, we assume / = f e is a combustion- type nonlinearity, f(u) = 
if u G [0,9]. Let e > be given, and consider the family of regularized 
problems 

' V° = eU: s + (kd s + V y ) 2 U a + ab- (kd s + V y ) U° + c a U° + f(U°) , 
< U a (s, y, r) 1-periodic in x and r, 

lim U a (s, y, r) = , and lim U a (s, y, r) = 1 , 

(2.1) 

where the parameter a G [0,1]. In this section, we will usually suppress the 
dependence of U a = U a,e,e and c CT = c a>ei g on e and on 9 since they remain 
fixed throughout the section. 

Before discussing the existence of solutions of (|2.1|) . let us point out some 
of the properties that solutions of (|2.1|) must possess. Since the equation in 
(12. lj) is nondegenerate, one may directly use the classical parabolic maximum 
principle to prove the following (cf. [2*§].|21j): 

• < U a < 1 on R x Q x (0,1). 

• U a (s,y,r) > U a (s',y,T) for s < s' and (y,r) GQx (0,1) (monotonic- 

ity). 

• The speed c CT > for which (j2.1|) has a solution {U a , c CT ) must be unique. 

• The solution U u must be unique modulo translation in s. 

• The solution U a converges exponentially fast to the limits as s — > ±oo. 

Our goal is to find a classical solution for this regularized problem when 
(7=1. It is known ([29 ,[21j) that there exists a solution of (|2.1j) when cr = 
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(i.e., in the absence of advection), so we apply a continuation argument 
to perturb this solution and solve (|2.1|) when a = 1. We now outline the 
application of this method to (j2.1|) ; the argument generalizes the arguments 
of [HI and |2U- 

Assuming that (j2.1j) has a classical solution (U°, cq) when o = ao, we seek 
a solution (U a ,c a ) of (|2.1|) for small perturbations of a > a . Substituting 
the expressions 

C(T = Co + ac a and U a = U° + aV a 

into equation (|2.1|) and using the fact that 

U° T = eU Q ss + (kd s + Vy) 2 U° + cx b ■ (kd s + V y ) U° + c U° s + f(U°) , 

we obtain the following nonlinear fixed-point problem for the perturbations 
Or and V a : 

V? - eVZ - (kd s + V y f V- - c V: - f(U°)V° = F(V% a, c a ) , (2.2) 

where the terms not appearing explicitly on the left have been absorbed into 
the nonlinearity F. 

To establish the existence of a solution V a of (J2.2j) . we employ the argu- 
ments used in [59j . Consider the weighted norm 

:= / u 2 pdsdydr, (2.3) 

JMxQx(0,l) 

where the weight p : R — > (0, oo) is defined by 

p(s) := 1 + e~ 25s 

for an appropriate constant < 8 < 1. The weighted spaces L 2 p := L 2 p (M. x 
Q x (0, 1)) and H l p := H p (R x Q x (0, 1)) are defined in the natural way. 
Define the simple function A : R — > R by 



f(l), if s>0, 
0, if s < 0, 

and let L and Li be operators with the common domain 

D(L) = D(L0 :={f 6 H p (EL x Q x (0, 1)) : (£;d s + V,) 2 «6^} 
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and respective definitions 

Lv := v T — ev ss — (kd s + V y ) 2 v — c v s — f'(U°)v and (2.4) 

Liv := v T - ev ss - (kd s + V y ) 2 v - c v s - A(s)v , (2.5) 

for v G D(L). Observe that L is the linear operator defined by the left-hand 
side of equation (|2.2|) and that we want to solve LV a = F(V a , a, c a ). 

As in [2T| and [2*9] . we now consider the invertibility of the operator L. 
Because of the regularization with e > 0, we can use the maximum principle 
and standard parabolic estimates to show the following: L\ is invertible on 
L 2 p , L\ and L differ by a relatively compact operator, and L is a Fredholm 
operator of index zero, whose adjoint L* has a 1-dimensional null space and 
an isolated eigenvalue of finite multiplicity at 0. Moreover, the eigenfunction 
v* of L* corresponding to is strictly positive, and the following Fredholm 
alternative holds: 

Proposition 2.2. Let f 6 L 2 p . The equation Lv = f has a solution if and 
only if 

/ fv* pdsdydr = 0. (2.6) 

</RxQx(0,l) 

When the solvability condition \2. 6]) holds, the solutions of Lv = f form a 
one-dimensional space. 

This result is a straightforward generalization of the corresponding results 
in [21] and (29J, so we do not include the details. Now, to solve equation 
(|2.2j) . we can use Proposition 12.21 to set up an iteration scheme, as in Section 
2.3 of "2Z], to which the contraction mapping principle applies. With each 
iteration, the constants {c n } < ^ =l that converge to c a are determined by the 
solvability condition ()2.6|) . As a result, we find that (j2.1j) is solvable for small 
perturbations of a: 

Theorem 2.3. If problem \2. 1)) has a classical solution (U°, Co) when o = gq, 
then there exists a positive constant a* = a*(U°,co) such that \2. 1)) has a 
unique classical solution (U a ,c a ) whenever a < a < a* . 

For details of this iteration procedure, see [21] and [2*7] . 

Next, having shown that solutions (U a ,c,j) of ()2.1|) exist for a < a*, we 
consider the behavior of solutions as a — > a*. We first prove bounds on the 
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sequence of speeds {c a } so that we can extract a convergent subsequence. 
Using periodicity and the conditions at infinity, integrating equation (|2.1|) 
over E x Q x (0, 1) yields 

c a = J f(U e7 )dsdydr>0, (2.7) 

from which we see that all of the speeds Cq- circ positive. Furthermore, as in 
[TUj . multiplying equation (j2.1j) by U u and integrating over IR x Q x (0, 1) 
yields the L? estimate 

J \W S +V y U a \ 2 dsdydr = -y + J U a f(U a ) ds dy dr < y , (2.8) 

where we used the identity ()2.7|) and the fact that < U a < 1 to bound the 
middle term. In addition, we have the following essential estimates on the 
speeds. 

Proposition 2.4. There exist constants c and c, independent of both a and 
e, such that 

< c < c CTi£ < c for all a E [0, a*), e G (0, 1] . (2.9) 

Proof. If ( V, c e ) solves the one- dimensional problem 

(1 + e)V ss + (a + c e ) V s + f(V) = , (2.10) 

with the appropriate limits at s = ±oo and the constant coefficient a defined 
by 

a := inf cr*6(u, r) • k < , 
Qx[o,i] vy ' ; - 

then the proof of Lemma 3.1 in |2H] shows that 

Ca,e < c e , for all a G [0, a*). 

Note that since a a constant, the constant 

a + c £ 



+ 



f2.11^ 



is the unique speed associated with the one-dimensional traveling wave prob- 
lem 

V ss + ~CeV B + rr^fiV) = . (2.12) 
(1 + e) 

9 



(See |2] for proof of existence and uniqueness of c e .) Because e G (0, 1], a 
comparison argument as in [2E] shows that c € is bounded above by c > 0, 
the unique speed determined by (j2.12j) with e = 0. Therefore, from (|2.11|) 
we have 

c e < (1 + e)c -a< 25 + trl&Hco = c, Ve G (0, 1]. 
This proves the upper bound in ([2.9)1 . As in [ID) , define the constant c by 



c := 



o 



We have 



Jy/2f(g)d£j = (^J^JW)k-(kU: + V y UndsdydT 

J f(U a )dsdyd^j ^2 J \W a + V y U° '\ 2 ds dy dr 



< 



< (O 2 , 

from which the lower bound follows. Here we have used ([2.8)1 and the fact 
that VI < for all (s, y, r) G R x Q x [0, 1]. 

□ 

Proposition 12 . 41 proves that the wave speeds are bounded away from both 
and +oo and thus have a positive limit point as a —>■ a* . Using this fact 
and standard parabolic regularity results, we obtain the following: 

Theorem 2.5. Let f be a combustion nonlinearity. For each e > 0, there 
exists a classical solution (U £ , c £ ) of 

U e T = eU e ss + (kd s + W y f U e + b- (kd s + V„) U £ + c £ U e s + f(U £ ) , 

U e (s,y,r) 1 -periodic in y and r, 

lim U £ (s, y, r) = , and lim U £ (s, y, r) = 1 . 

s— >oo s— >— oo 

(2.13) 
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The solution U 6 is unique modulo translation in s. Moreover, 

< U £ < 1 on R x Q x [0, 1] , 

U £ (s,y,r) > U £ (s',y,r) for s < s', y G Q, t 6 [0, 1] , and 

c £ > 0. 

Proof. Fix an e > 0. By interior parabolic regularity results (^S]) and the 
bounds on / and U a , 

l|£> 2 £/1U + ||t/?IU < c (II^IU + < C (1 + ll/IU) |£| , 

(2.14) 

where CC KxQx (0,1) and the constant C is independent of a. As 
the upper bound on the right is independent of a, choosing p sufficiently 
large ensures that U a is Holder continuous. Schauder estimates then apply, 
and the Arzela-Ascoli Theorem guarantees that U a — > U* uniformly on any 
compact set as a — » a*, where U* solves the equation in (|2.1|l for cr* (here 
we pass to a convergent subsequence, if necessary). Since problem ()2.1|) has 
a solution at cr*, solvability may actually be continued all the way to a = 1, 
yielding the solution U e of problem ()2.13|) . Because the convergence of the 
functions U a is uniform, the boundary conditions at s — ±oo may be verified 
as in [21] and and the monotonicity and uniqueness properties follow 
from the classical parabolic maximum principle, as described above. 

□ 

2.2 Regularized Equation: Positive Nonlinearity 

Now that we have existence of solutions to the regularized problem (|2.13J) 
for the combustion nonlinearity we wish to show existence of solutions to 
()2.13|) for a KPP nonlinearity /. As in [3J, we approximate / by a sequence 
of combustion nonlinearities as follows. Given 9 G (0, 1), let \ 9 '■ K — > [0, 1] 
be a smooth cut-off function such that 

' x \u) = 0, for u < 9/2, 

< x 6 (u) = 1, for u > 9, and (2.15) 
X e >x e ' if 0<9'. 
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For each 6, the function / := \ f is a nonlinearity of combustion type, and 
the increasing sequence f 9 converges uniformly to / on [0, 1]. 

For a fixed e > 0, the results in the preceding section show that, for each 
9, there is a unique classical solution (U 9 , eg) of 

' U 9 = sU° ss + (kd s + V,) 2 U e + b ■ (kd s + V,) U 9 + c e U° s + f 9 (U e ) , 

< U (s, y, t) 1-periodic in y and r, 

lim U (s, y, r) = , and lim ?7 (s, y, r) = 1 , 

S— >00 S^r — OO 

(2.16) 

which satisfies the normalization condition 

min{f/ e (0, 2 /,r)} = i (2.17) 

y,r 2 

Using the parabolic maximum principle, one can show that the wave speeds 
eg increase as 9 approaches 0, since the sequence f 9 increases as 9 decreases 
( 29j). Moreover, these speeds satisfy 

eg < c = eg , for all 9 > , 

where eg was defined above in the proof of Proposition 12.41 Although the 
constant eg depends on 9, the estimates of Section 8 in [3j show that the 
speed of the corresponding one-dimensional traveling waves 1)2.12)1 can be 
bounded above independently of 9. Hence, the constants eg are bounded 
independently of 9. The sequence eg therefore converges to c* > as 9 — >• 0. 

Having shown that the speeds eg converge as 9 — > 0, we can apply 
parabolic regularity (as in the proof of Theorem 12.5)) to show that the func- 
tions U 9 converge, as 9 — >• 0, to the unique classical solution U £ of the problem 

' U £ T = eU e ss + {kd s + V y f U £ + b- {kd s + V„) U £ + c* e U £ s + f{U £ ) , 

< U £ (s,y,r) 1-periodic in y and r, 

lim U £ (s, y, r) = , and lim U £ (s, y, r) = 1 , 

s— >oo s^— oo 

(2.18) 

subject to the normalization (j2.17j) . The fact that the limits hold at s — ±oo 
follows from the normalization ()2.17|) and arguments as in [21] that make use 
of local uniform convergence. 
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2.3 Removal of the Regularization 

The previous section established that, for any given e > 0, there is a classical 
solution (U £ , c £ ) of problem ()2.13|) for the positive nonlinearity. To obtain 
a solution of the original problem ([1.7)1 . we want to let the regularization 
parameter e go to 0. Therefore, we must prove that the speeds c £ and the 
functions U £ converge in some sense as e — » 0. Proposition 12.41 shows that 
the speeds c £ are bounded away from and +oo, while the next result yields 
the needed estimates on the functions U £ . 

Proposition 2.6. The solutions U £ of are uniformly bounded with 

respect to e in W l( ^. 

Proof. Since < U £ < 1 for all e > 0, U £ G V\ oc for any p > 1. It follows 
from the monotonicity of U £ with respect to s that U £ s < for all e, and we 
have 

U e s dsdydr = -\ for all e, (2.19) 

5xQx(0,l) 

showing that U\ G L\ oc . Multiplying the equation in ()2.13|) by U £ and inte- 
grating (as in the derivation of ()2.8|) above), we have 

J \kU £ s +WyU e \ 2 dsdydr = ~ + J ' U £ f(U £ ) ds dy dr < ^ < C - , (2.20) 

for the constant c from Proposition 12.41 It follows from ()2.2()|1 that 
\kU e s +W y U £ \ G Lj oc , and we conclude from (JQH) that \W y U £ \ G Lj oc . Fi- 
nally, rearranging equation ()2.13|) . multiplying by (U £ — c e U £ ), and integrat- 
ing, we have 

\U £ T - c £ U £ \ 2 dsdydr = J b ■ (kU £ + V y U £ ) (U £ - c £ U £ ) dsdydr . 

Applying Holder's inequality and using ()2.8j) . we see that (U £ — c £ U £ ) belongs 
to L 2 , and it follows from (J2.19j) that U £ G L\ oc . We therefore find that 
U £ G and it is clear that, for any compact E, the W 1,1 (E) norm of U £ 

can be bounded by a constant independent of e > 0. 

□ 

Since the functions U £ are bounded in W^ c , there exists a limit function 
U G L\ oc (in fact, U G -BV[ oc [Oj) such that, after passing to a subsequence, 

U £ — >U in L\ oc and U £ — >U a.e. . 
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Of course, the U £ also converge weakly in Ll oc to U ; letting e — > in equation 
(|2.13|) . we see that U solves the traveling front equation (|1.8|) in the sense of 
distributions. 

To analyze this convergence more closely, define the functions u e by 

u £ (x, t) := U e (k ■ x — c £ t, x, t) . 

The uniform L 2 estimates on \kU £ + V y U e \ and \U £ — c £ U £ \ show that the 
sequence u £ is compact in L\ oc and thus converges (up to subsequence) to 
u G Hl oc . It follows that a is a weak solution of the original semilinear 
equation (jl.lj) . and parabolic regularity results then imply that u is a classical 
solution of p. 1)1 . The directional derivatives 

U T - cU s , k% + U yi , % = 1, . . . , N, and (kd s + V ' y f U 

are therefore smooth, but we cannot conclude that U is a classical solution 
of the traveling front equation (jl.8j) . 

Now we show that the function U satisfies the desired limits s — > ±oo: 

]jmU(s,y,r) = Q and lim U(s, y, r) = 1. (2.21) 

S— HX S—f — OO 

The above bounds on U € are independent of translation in the s direction, 
so we are free to normalize the sequence U e by 

/ U e {s,y,r)dsdydr =-, for all e > 0. (2.22) 

7(0,1] xQx [0,1] 2 

Since the sequence U e converges to U in Lj oc , this normalization condition 
also holds for the limit U. Since U\ < 0, we conclude that U is nonincreasing 
in s: 

U(s,y,r)>U(s , ,y,r) iis<s', V(y,r) G Q x [0, 1]. 
Therefore, since we also have < U < 1, the functions 

^ + (y,T)= lim U(s,y,r) and if)~(y,r)= lim U(s,y,r) 

are well-defined. We claim that ip + (y,r) = and ip~{y,T) = 1. To see 
this, let h(y,r) be any smooth function periodic in (y,r) G Q x [0, 1]. For 
n = 1,2,3..., define the sequence 4> n (s,y,r) = ^ n (s)h(y,r) where £ n (s) = 
^o(s — w), and £ satisfies 

< Co(s) < 1, for all s G -R, 

^o(a) = 0, for |s| > 1, (2.23) 
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and 

/ £ Q (s) ds = I. 

JR 

Note that for each n, £ n satisfies 

(Q s ds = 0= [ {Q ss ds. (2.24) 

l JR 

Now multiply ()1.8|) by and integrate over Rx Q x [0, T]. We see that 
/ U(s,y,r)(4> n ) ss dsdydT = h(y,r) ( U(4> n ) ss ds) dydr 

JRxQx[0,l] JQx[0,l] \Jr J 



h(y,r)( / ty + + (U-ij> + )](<t> n )„ds) dydr. 
x[0,l] \Jr J 

For fixed (y, r), we have 

lim \U(s,y,T)-iP + (s,y,T)\ = 0. (2.25) 

s— >+oo 

Therefore, for each (y, r), we have 

lim f [U -i) + }{(j) n ) ss ds = 0. (2.26) 

n->oo 

Also, by definition of (j), for each (y, r), 



^ + (y, r)(0") ss = V + (2/, r) / (0") ss rfs = 0. (2.27) 

R JR 

Using the dominated convergence theorem, we now conclude that 

lim / U(<p n ) ss dsdydT = 0. (2.28) 

JrxQx[0,1] 

Applying a similar argument to each term in the equation, we let n — > oo 
and conclude that ip + satisfies 

J tp + h T dy dr — J ip + A y hdydr — J tp + b ■ V y h dy dr = J f(i/) + )hdy dr. 

(2.29) 
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Since h was arbitrary, we conclude that ip + is a distributional solution of the 
periodic-parabolic equation 

Vv - A y ip + b ■ V„V = /(V) (2.30) 

Because of the bounds on [/, we know that < ip + < 1, so that right hand 
side f(ip + ) G L P (Q x [0, 1]), for any p > 0. By parabolic regularity estimates, 
it now follows that tp + e H 72 ^ for any p > 0. Then, as in the analysis 
following ([2. 14)1 . we find that ip + is a classical solution of ([2.30)1 . assuming 
b is sufficiently smooth. The maximum principle implies that ip + must be 
constant. Hence, f{ip + ) = 0, so that tp + = or ip + = 1. Because of the 
normalization ()2.22|) and because U is nonincreasing, ip + < ~. Therefore, 
ip + = 0. A similary argument with £ n (s) = Co(s + n) implies that ip~ = 1. 

We have now shown that U satisfies the desired limits, concluding the 
proof of Theorem 12.11 

3 Propagation Speed for KPP Nonlinearity 

In this section, we make the additional assumption that / is a KPP non- 
linearity satisfying f(u) < uf'(0). Under this assumption, we find that the 
speed c* of the traveling wave obtained in the preceding construction is min- 
imal and can be characterized by a variational formula which extends the 
well-known variational formula for the KPP minimal speed when b is time- 
independent. Moreover, we can describe propagation of front-like solutions 
at any speed higher than c*. First, we characterize the speed c*. 

3.1 Characterization of the KPP Minimal Speed c* 

To introduce the variational characterization, we define //(A) to be the prin- 
cipal eigenvalue of the periodic-parabolic operator ([T2~]) 

L A $ := A x .$ + (b - 2Xk) ■ V x $ + (A 2 - Xb ■ k + /'(0)) $ - $ t . (3.1) 

defined on spatially-temporally periodic functions $(x,t). The associated 
eigenfunction $ = $a > is unique up to multiplication by a constant. We 
have the following characterization of the speed c* obtained in the preceding 
construction: 
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Theorem 3.1. The speed c* of the traveling wave obtained in the preceding 
sections satisfies 

e = inf ^ (3.2) 

A>0 A 

Moreover, this speed is minimal in the sense that if (U, c) is any traveling 
front solution of 11. 8\) . we must have c> c* . 

Proof. In the special case that b is a time-dependent shear flow, this char- 
acterization was proven in [25- The proof there generalized ideas of [2j and 
relied heavily on regularity properties of the traveling wave equation. In the 
present case, the traveling wave equation is degenerate. Nevertheless, for the 
regularized equation, the arguments of [2Q can be extended in a straightfor- 
ward way to prove that 

-S-gf^, (3.3) 

where n £ (X) denotes the principal eigenvalue of the periodic-parabolic oper- 
ator 

L £ ' A $ := A x <& + (b- 2Xk) ■ V a .$ + ((1 + e)A 2 - Xb ■ k + f(0)) $ - $ t , (3.4) 
defined on spatially-temporally periodic functions t). Note that 

//(A) = fi(X) + e\ 2 . (3.5) 

so that 

iBi m =iBf (m +eX ) (3. 6) 

a>o A a>o y A / 

Using the properties of /i(A) described in the appendix, we see that the speeds 
c* converge to 

c* = inf ^ (3-7) 

A>0 A 

as e — >• 0. We will show below that this limit c* is indeed the minimal speed 
at which a traveling wave solution of (jl.lj) may propagate. □ 

3.2 Generalized Propagation and Minimality of c* 

In this section we show that the speed c* defined above is indeed the minimal 
speed of propagation and that solutions of the corresponding initial value 
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problem can propagate at speeds higher than c*. To do so, we generalize 
some results found in [Tj and jTHl that define the asymptotic propagation 
speed in the KPP case, without relying on the existence of smooth traveling 
waves. We consider propagation in a given direction k, a unit vector in R N , 
arising from solutions to the initial value problem 

Lu + f(u) = A x u -u t + b(x, t)-Vu + f(u) = (3.8) 

u(x, 0) = Uo(x). 

where / is the KPP nonlinearity satisfying f(u) < f'{0)u for u G (0, 1). For 
a given r G R, we use S~ C R N to denote the open half-space: 

T,~ = {x E R N \ k-x<r}. 

Similarly we define by (k ■ x) > r, i.e. the interior of the complement of 
E~. Also, we define for r > 0, 

S° = {x E R N | \k-x\ < r} . 

Unless otherwise stated, we will consider initial data uq satisfying < no < 1 
and the limits 

lim inf u (x) = 1 and lim supM (x) = (3.9) 

That is, Mo will converge to limits and 1 uniformly in the directions k 
and —k, respectively. In this case we say Uq is wave-like in the direction 
k. The main tool for this analysis is the eigenvalue problem associated to 
propagation in the direction k: 

L\(j) = A(f)-(f>t + (b- 2\k) ■ V0 + (A 2 - A(6 ■ k) + /'(O))0 = n(\)<j> (3.10) 

with <ft(x, t) = <p\ > being periodic in both x and t. This is obtained by 
substituting a solution of the form e'^^'^cp^x, t) into the linearized version 
of equation ()3.8|) and taking /x(A) = Ac. Define c* by 

c * = iri f^ >0 . (3.11) 

A>0 A 

Let A* > be the point where this infimum is attained. As described in 
the appendix, for each c > c* there exists a unique A c G (0, A*) such that 
/i(A c ) = cA c . 

Note that c* and A* depend on the direction vector k. If we consider 
propagation in the opposite direction, k — > —k, then (|3.11|) gives us another 
c** and A** corresponding to —k, and it is not necessarily true that c** = c*. 
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Definition 3.2 (Speed of Propagation). We say that the solution to the 
initial value problem \3. 8\) propagates in direction k with speed c if for any 
d > c and any r G R, 



lim sup u(x + c'tk, t) = 0, 

and for any c" < c and any r G R, 

lim inf u(x + c"tk, t) = 1. 

The following theorems and their corollaries allow us to characterize the 
asymptotic speed in terms of the decay rate of the initial data. Here are the 
main results: 

Theorem 3.3. Let u satisfy \3. 9i) and u (x) < Ce~ Xck ' x for some constant 
and c> c* . Then for any d > c and any r e R we have 

lim sup u(x + dtk,t) = 0. (3-12) 
That is, the solution u propogates with speed < c. 

Theorem 3.4. If there exists a constant m > and an interval [a±, 02] C R 
such that 

< u (x) < 1 for all x G R N , 
u (x) > m > wherever (k ■ x) G [01,02], (3.13) 



then 



lim inf u(x + c'kt,t) = 1 (3.14) 



for any d G (— c**, c*) and r G R. 

Theorem 3.5. Let c> c* and let uq satisfy \3. ty) and uq(x) > C\e~ Xck x for 
[k ■ x) > C 2 , for some constants Ci,C 2 - Then for any d < c, we have 

lim inf u(x + c'tk,t) = 1. (3.15) 
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Corollary 3.6. For initial data satisfying \3. Mj) . 

lim inf u(x + c'tk, t) = 1 (3.16) 

for any c' < c* . 

Corollary 3.7. If uq satisfies \3. 9\) anduo(x) < Ce~ x * k ' x for some constant 
C > 0, then the asymptotic speed of propagation in direction k must be c*. 

Recall that — c** corresponds to propagation in the — k direction at speed 
| — c** | . The point is that an initial disturbance must spread at a rate at least 
as fast as the minimal speeds in the corresponding direction. In what follows, 
we rely heavily on the eigenvalue problem (j3.10j) to construct global sub- and 
super-solutions to the initial value problem (|3.8j) . 

Proof of Theorem 13.31 For any c\ G (c, c'), we have C\\ c > /x(A c ), by 
definition of A c . Then define ip(x,t) = Me~ x °( k ' x ~ Clt '(f)(x,t), where = 4>\ c > 
is defined by (|3.10p . and the constant M is chosen so that Wo(^) < ip{%, 0) 
for all x G R N . Then ip is a supersolution: 

Lip + f(0)ip < 
ip(x, 0) > u (x). 

Since f{ip) < tpf'(0), it follows from the maximum principle that ijj(x,t) > 
u(x,t) for all x G R N , t>0. Therefore, 



lim sup u(x + c'tk,t) < hm sup ip(x + c'tk,t) 

= lim sup Me- Ac(fc -* +(c '- cl) *V(z + c'tk, t) 
< M||0|| oo e- Acr lim e - Ac(c '- Cl) * = 0, (3.17) 



t^QO 



since d > c\ . This proves Theorem 13.31 

The proof of Theorem 13.41 extends ideas from Theorem 3.4 of First, 
we state the following lemma, which is a straightforward extension of Lemma 
3.5 in P2J: 

Lemma 3.8. There exists a d < c* such that for all c G (c',c*), there exists 
A eC\TZ and a function <p = <f>(x,t) G C 2 > l {T n+l ,C) with Re((j)) > and 
solving 

A</> - fa + (b - 2Xk) ■ V(f) + (A 2 - Ac - A(6 ■ k) + f'(0))<f> = 0. (3.18) 
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The important point of the lemma is that we can choose A to be strictly 
complex and to have strictly positive real part. This will allow us to 
construct global subsolutions that have no tail (or sufficiently small tail) and 
can fit beneath uq(x). We omit the proof of this lemma since it is identical 
to that of Lemma 3.5 of JH], except that the operator is somewhat different 
because of the time-dependence and more general flow structure. 
Proof of Theorem 13.41 

Step 1: If we replace f (0) by f(0) - 5 in (jB~TUJ) . then from (PHTI) we 
obtain a new minimal speed c* 5 instead of c* with c| /* c* as 5 — > 0. Also the 
lemma would hold with /'(0) replaced by /'(0) — 5, so that for S sufficiently 
small and for c R < c* 5 < c* sufficiently close to c|, we apply the lemma to 
obtain a function such that 

A0 - <f> t + (b - 2Xk) ■ V0 + (A 2 - Xc R - X(b ■ k) + f (0) - 5)0 = 0, 

where Re((p) > and A G C \ 1Z depends on cr. Now let 

i/j = Re{e- x{k - x - CRt) (f){x,t)). 

Then if) solves 

^ _ - b ■ W = (/'(0) - 5)^, 
or in the moving frame ip(x, t) = ip(x + CrH, t): 

i>t - - (b + c R k) ■ = (/'(0) - 5)^. 

Now 

ip(x, t) = e - Xr{k - x - CRt) [<j) r cos(\i(k ■ x - c R t)) + <f>i sin(Ai(A; • x - c R t))\ 

^ (3-19) 

where A r , Aj, <f> r , (pi denote the real and imaginary parts of A and 0. If (k ■ x — 
c R t) = ^ for n G Z, then ip > 0, since r > 0. Similarly, if (k • x — c R t) = 
(2n+i)n ^ faen ip < 0. It follows from i.e. from the complexity of A, 

that for a given rj > 0, there exists an unbounded set D C R N+l and an 
interval (ai,a 2 ) C -R such that 

< ip(x, t) <rj for (x, t) G D, 
ip{x,t) = for (x,t) G 
C {(x, t) | {k ■ x — c R t) G (ai, a2)} • 
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Furthermore, we may choose some a 3 G (a\, a 2 ) such that a 3 = ^ f° r some 
n E Z. Now we define 

if G D 

(3.20) 

otherwise. 

So, < ij)R < 7] with support D bounded by the two hyperplanes (k-x—cut) = 
a\ and (k • x — cr£) = a^. Also, i[)r(x, t) > whenever (k ■ x — CRt) = a 3 . If 
we choose r\ > such that f(s) > s(/'(0) — 5) for s G (0, rj], then satisfies 

(Vfl)t-AVji-&. </(V>h) 

for all (x, t) G -D. By translating ^ in the variable t, we may assume that 
has support in the (ai,a 2 ) defined in the statement of the theorem, so 

that by reducing r], we may assume iPr(x, 0) < u (x) for all x G R N . 
Now ipR is the desired subsolution. The function u(x,t) satisfies 

u t — Au — b • Vu = f(u) 
u(x, 0) = u (x) > i/;r(x, 0). 

By maximum principle u > for t > 0, and w > for all t > 0. 

Suppose we repeat the preceeding analysis after switching k — > — /c. In 
this case, we have a new minimal speed c|* corresponding to waves in the 
direction —k. Furthermore, we can create a bump function ipi (analogous 
to %1)r) that is a subsolution propagating in the — k direction with speed cl 
arbitrarily close to c* 5 * and — c}* < —cl < cr < c}. It follows that if u(x,t) 
satisfies 

u t — Au — b ■ Vw = f(u) 
u (x) > i/Jr(x, 0) and u (x) > ^ L (x, 0) , 

then u > ipR and u > ipi for all t > 0. In particular, there exists e > 0, A > 
and constants a° L , a° R such that 

u(x, t) > ijj L {x,t) > e > if \k ■ x — a L (t)\ < A 
u(x,t) > ip R (x,t) > e > if \k-x- a R (t)\ < A, (3.21) 

where a L (t) = a° L — c^t and a R (t) = a° R + crL The constant e is determined 
by r] and inf^,*) (fi r ,R, and in.f( x ,t) (fir,L, where <fi r ^R and (p r< L are the real parts 
of defined in the construction of ipR and ipL, as in (j3.19j) . 
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Step 2: Now with az,(t), or(£) defined as above, we define for r > the 
set n T C R N+1 by 

Q T = {(x,t)\ k-xe [a L (t),a R (t)],t> r} . 

and we show that 

lim inf«(a;,t) > e. (3.22) 

In other words, u stays bounded uniformly away from zero in a region spread- 
ing at rates cl and c R in the directions — k and k, respectively. 

Let / G C 1 ([0, 1]) be a nonlinear function satisfying f(u) < f(u) for all 
u G [0, 1] and 

f(u) = for uG {0}U [e, 1], 
f(u) > for u E (0,e). 

Now let to > such that a^i) < a^(t) for t >t and take 
g = -inf {u(x,t )\ (k ■ x) e [a L (t ),a R (t )]} 

Z x 

Let the function g(t) solve g'(t) = f{g(t)) for t >t and g (to) = go- Then in 
the region Q tQ , we have 

g t -Ag-b-Vg = f(g)<f(g) 
u t — Au — b ■ V« = f(u) 

with u(x, t) > g{t) in the parabolic boundary dQ to , since by (j3.21|) u > e on 
the lateral boundary of f2 to . Maximum principle now implies that u(x,t) > 
g(t) for all (x,t) G fi to . Since g(t) — > e, this implies the ()3.22j) . 

Step 3: In this step we improve the result of Step 2 by showing that 
for any q G (0, 1/2) arbitrarily small, then for some constant d = d(q) to be 
determined, 

lim mfu(x,t) > 1 - q (3.23) 

r^oo no- 
where 

Q' T = {{x,t)\ k-xe \a L {t) + d,a R (t) -d],t>r}. 

Inspired by the technique of [7] (see p. 158), we let 9 G (0, 1) and define the 
linear function 

f(u) = s(i-e-u) 
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choosing the constant s > so that f(u) > f(u) whenever u > e. Thus by 
the preceeding step, there is t x > sufficiently large so that 

u t — Au — b ■ Vu > f (u) for all (x, t) G Vt tl . 

For constants 7, k G (0, 1) to be determined, define 

1 

a L (t) = a L (t) + ( ^ 4 
a R (t) = a R {t) - ( - 
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and the function 



k, if (k ■ x) G [a L (£),a/?(t)] 
p(x,t) = <{ -7(fc • x - a R (t)) 4 + «, if (A; • x) > a fl (f) (3.24) 
— 7 (A; • x — a L (t)) 4 + k, if (A; • x) < a^t) 

Note that the functions a L {t) and a R (t) were defined so that p(x,t) > if 
(x,t) G £l t and p < otherwise. Choose er G (0, s) and let 

g(t) = 1-9- e- a{t - h) 

Now we compare with u the function defined by 

z(x,t) = g(t)p(x,t) 

We compute: 

z t - Az-b-Vz- f(z) 
= g\t)p + g(t)p t - g(t)(Ap + b ■ V p) - s{l - 6 - gp) 
= (a- s)e-^-^p + g(t)(jH -Ap-b-Vp) + s(l - 6){p - 1) 
< (l-9)(p t - Ap-b- Vp) + s(l-0)(/c-l) (3.25) 

A straightforward computation shows that if (x,t) G Qt, where p > 0, then 
there is a constant C > such that 

|/»t| < |Vp|<C 7 3, \Ap\<C^ 
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Thus choosing 7 sufficiently small, we can make (1 — 9)(p t — Ap — b • Vp) 
arbitrarily small, so that z t — Az — b ■ — f(z) < since k < 1. Therefore, 
in the region Q tl , 

z t - Az-b-Vz< f(z) 

ut - Au - b ■ Vu > f(u) (3.26) 

with u(x,t) > z(x,t) on the parabolic boundary dQ tl . Maximum principle 
implies that u(x,t) > z(x,t) for all (x,t) G Q tl . Therefore, 

lim inf u(x, t) > lim inf z(x,t) = k(1 — 6) (3.27) 

where 

tf T = l{x,t)\ k-xe[a L {t)+ ,a R {t)- (^y\,t>r\. (3.28) 

Since 9 and k can be chosen arbitrarily close to and 1, respectively, we can 
make «(1 — 9) > 1 — q, proving the claim ()3.23|) . 

Step 4: Now we complete the proof. Recall that a^(t) = a° L — c^t and 
a n(t) — a R + c Rt- Therefore, the result ()3.23|) of preceding step implies that 
for any d G (— Cx, cr) we have 

lim inf u(x + c'kt, t) = 1 (3.29) 

By definition of cl and cr, we can choose cl and cj? arbitrarily close to c* 5 * 
and c|, respectively. Furthermore, we can choose 5 sufficiently small so that 
Cg* and c* s are arbitrarily close to c** and c*, respectively. Hence the result 
(I3.29|) applies to any c' G (— c**,c*). This completes the proof of Theorem 

EH 

Proof of Theorem 13.51 Without loss of generality, assume C2 = 0. Thus, 
uo{x) > Ce~ Xck ' x for (k ■ x) > 0. We now construct a subsolution of the form 

i/)(x,t) = d ie - Xlik - x - C2t) <p Xl {x,t) - d 2 e- X2{k - x - C2t) <p y (x,t) (3.30) 

for some Ai,A2 satisfying A c < Ai < A2. As described in the appendix, if 
we replace f'{0) with f'(0) — 5 in the eigenvalue problem, minimal speed c* 
becomes c}, with c* 5 < c*, c* s — > c* as 5 — > 0. Let ps(\) denonte the eigenvalue 
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corresponding to the problem with /'(0) — 5. Then, we can pick 5 small so 
that there exists Ai, A 2 and c 1; c 2 such that 

MAj /i^(Ai) 

c 2 = — < — = ci (3.31) 

A 2 Ai 

and c 2 < ci < c and A c < Ai < A 2 (see appendix for properties of fig). Now 
define the function 

<ijj(x, t) = d ie - Al(fc - x - C2t Vi(2;, t) - d 2 e- X2(k - x - C2t) (f) 2 {x, t) (3.32) 

where 0i,0 2 > are the eigenfunctions defined by (j3.1()j) with /'(0) — 5, 
corresponding to Ai,A 2 . Given an r\ > 0, we can pick d\,d 2 > such that 
ip(x, 0) < Uq(x) for all x G -R^ and sup^^ ift(x,t) < rj. This is possible since 
ti (x) > Ce~ Acfc ' x for (fc-x) > and u satisfies f|3.9|) and A c < Ai < A 2 . Also, 
for r sufficiently large, 

?/>(x, 0) > for x G £+. 

By (|3.31|) . Aic 2 0i < fj,s(^i)(j>i- As a result, ^ is a subsolution to the linearized 
problem: 

^-A^-b- VV< (/'(0)-o> 
V>(x,0) < w (^)- 

Now, since /(s) > s(/'(0) — 5) for s G (0, r/) and since t/> < for all (x,t), 

Vt-AV-6-V^</(V) 

?/>(x, 0) < Wot^)- 

By maximum principle, u > ip for all (x,t). As in (|3.21|) . we find conclude 
that there is an e > and A > and a constant ao such that 

u(x,t) > 4>(x,t) > e > if \k ■ x - a R (t)\ < A 

where or(£) = ao + c 2 t. 

Now using the assumption that 

lim inf Uq{x) = 1, 
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we proceed as in Step 2 of Theorem 13.41 to show that 

lim inf u(x, t) > e 
t — >oo ri T 

on the set 

tt T = {{x,t)\ k ■ x E (-00, a R {t)],t > t}. (3.33) 

Then it is easy to extend the arguments of Step 3 and Step 4 in the proof of 
Theorem 13.41 to conclude the proof that 

lim inf u(x + c'tk, t) = 1. (3.34) 

Proof of Corollary 13.61 Suppose the initial data is wave-like, i.e. satisfies 
()3.9|) . Without loss of generality, we may assume that uq{x) > if (k-x) < 0. 
Then let v (s) > be a smooth function with compact support in the interval 
s G (0, 1), with \\v ||oo = 1. Then consider the function ipo(x) defined by 

00 

V'o( a; ) — /J a,jv(k ■ x + j) (3.35) 

where the coefficients {aj} are chosen such that ip(x) < Uq(x). Applying the 
argument from Theorem 13.41 the solution to the initial value problem 

Lij + fty) = 

V(x,0)=^ (3.36) 

must satisfy 

lim inf ijj(x + c'tk, t) = 1 (3.37) 

for any d < c*. By maximum principle, the result holds for u(x,t) since 
u {x) > ip . 

Proof of Corollary 13.71 This follows from the preceeding corollary and 
from Theorem 1.1, with A c = A*. 

4 Appendix: Properties of the Eigenvalue Prob- 
lem 

The eigenvalue problem associated with propagation in the direction k is: 
L\(f) = A0 - <j) t + (b - 2\k) ■ V0 + (A 2 - A(6 ■ k) + f(O))0 = /i(A)0 (4.38) 
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with (j){x, t) = (p\ > being periodic in both x and t. In this section we 
describe the following properties of the principal eigenvalue /i(A) as a function 
of A: 

Lemma 4.1. The map A i— > //(A) strictly convex in A. 

Proposition 4.2. There is a point A* G (0, +oo) at which the curve A i— > 
achieves a unique global minimum (over (0, +oo)j which is also the only 
local minimum. Hence, is strictly decreasing over interval (0, A*) and 
increasing over interval (A*, +oo). 

Proof of Proposition 14.21 Integrate the equation (jH.lOj) over D = Q x 
[0, T] and divide by j D <fi dxdt > 0. Periodicity and the fact that V ■ b = 
imply that 

MA) = A 2 + f(O) + A /D( r 6 '! )0 - ( 4 - 39 ) 

Since > 0, the last term is bounded by ||6||oo> so that /x(A) grows quadrat- 
ically as A — > ±oo. Letting A = 0, maximum principle implies that /i(0) = 
/'(0) > with eigenfunction being a positive constant. Then the proposition 
follows immediately from the fact that /i(A) is convex in A. 

The convexity of fi(\) follows from a minor extension of the argument 
in [2] (see Section 5.2 therein), but for convenience and since the eigenvalue 
problem is nonstandard, we include the details here. First, we make use of 
the following characterization of the principal eigenvalue, which holds even 
though the operator is not self-adjoint. Let operator L be defined by 

Lcj) = (i,,o r x , - 4>t + bi<p Xi + c0 (4.40) 

where = (aij(x,t)) is uniformly elliptic, hi = bi(x,t) and c = c(x,t), 
with summation notation implied. All coefficients are periodic with respect 
to both x and t. Denote by /i the principle eigenvalue of L over the space 
E = Cp(Q x [0, T]) of C 2 functions periodic in D = Q x [0, T}. 

Theorem 4.3. The following variational formula holds for the principal 
eigenvalue fi: 

\x = inf sup — — (4-41) 

( x>t )eD W 

where E + = {ip e E \ if> > 0}. 
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Proof of Theorem I4.3t This is an extension of the version proven in 
[2] for time-independent operators. Clearly when ip = <p is the eigenfunction 
for L, then = fi. Suppose there is a function tp > 0, ip E E such that 



< /i for all (x, t) G .D. 



(4.42) 



Then Lip — flip = m < for some function m. Now there exists r such that 
<p < rip for all (x,t) with equality holding at a point. This follows from 
positivity and continuity of ip and (p. Define the function w G E by w = ±. 
Plugging (p = wip into the equation Lcp — ficp = 0, we obtain 

«,j tr -w t + hw Xz + aijW Xi —^- + a ij w X] ^- = --r{^> ~ aWO > °- 

By definition of w and t, w < r with equality holding at one point. By 
maximum principle, we must have w = r. Note that periodicity in t is 
necessary to apply the maximum principle in this way. Hence <p> = rip, which 
contradicts ([4.42)1 . Hence 

sup — > fi (4.43) 

(x,t)eD W 

for all ip G E. This completes the proof. 

Proof of Lemma I4.lt The Lemma extends an argument in [2|, Section 
5.2; we include the details for convenience. From the theorem we have 



For a function ip G E + , we define ip = e ip. Then it is easy to check 
that 



V>££+ (x,t)€D V 

with the set E$ defined by = [e~ xk - x ip | ip e E + }. For A, 7 > 0, let 
(p G E~£ and ip G E+. If we define w = y/(jnp, then w G E£, where a = ^y 1 . 



fl(X) = inf sup 

i>&E + ( x ,t)£D W 



(4.44) 




(4.45) 
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Therefore, by (jQ5jl 



Let us compute 



M^PO < /'(0) + sup — . (4.46) 

^ (x,t)eD w 



Lw 1 Lcj) 1 



Note that 



1 /|V0| 2 |V^| 2 V0-WA V0-V^ 

4 V <P 2 V H ) H ' 

V0- W _ |V0| 2 _ |V?/;| 2 _ 
<p1p <p 2 1p 2 

(2(MW ■ W - V0| 2 - ^ 2 |V^| 2 ) . (4.48) 

If we take a = (p Xi ip and b = ip Xi <p and apply the inequality lab < a 2 + b 2 , 
we conclude that the expression in ()4.48|) is nonpositive. Now returning to 
(I4.47|) . we conclude that 



Therefore, 



Lw lL(p 1 Lip 
~w~ ~ 2~<f + 2~' 

Lw 1 Lcp 1 Lip 
sup < -sup — + -sup— . (4.49) 

{x,t)€D w 2 <p 2 ip 

Since <p and ip were arbitrarily chosen, we conclude from ()4.45j) and ()4.4(ij) 
that 

Therefore, A i— > fi(X) is convex. This completes the proof. 

Remark: It follows easily from ()4.45|) that if /xa(A) denotes principal 
eigenvalue of (IPK|) with f(O) replaced by f(0) - 5 for 5 > f'(0) > 0, we 
must have 

• /j, s (\) = yu(A) - 5 for all A > 

• ^j^- locally uniformly as 5 — > 

• also satisfies the properties described in Lemma f4. II and Proposi- 
tion 
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5 Conclusions 



We proved the existence of KPP traveling fronts propagating at minimal 
speeds in space-time periodic mean zero incompressble flows. The minimal 
speeds are given by a variational principle defined through principal eigen- 
value of a linear space-time periodic parabolic operator. Front existence holds 
also for general positive nonlinearity. It is interesting to study of the quali- 
tative properties of minimal speeds as the advection field varies its intensity 
in future works. 
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